Two-band model of Raman scattering on electron-doped high-T c superconductor 
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We have analyzed the B\ g and Big Raman spectra of electron-doped cuprate superconductors 
Nd2- a; Ce a :Cu04 and P^-nCe^CuCU using a weakly coupled two-band model. One of these two 
bands is centered around (±7r/2, ±7r/2) and couples more strongly with the Big mode, while the 
other is centered around (±tt, 0) and (0, ±tt) and couples more strongly with the B\ g mode. This 
model explains in a natural way why the Big Raman peak occurs at a higher frequency than 
the Big one at optimal doping, and how these two peaks change with doping in agreement with 
experiments. The result thus supports that there are two kinds of quasiparticles in electron-doped 
cuprates and d a .2_ H 2-wave superconductivity is driven by the holelike band and a proximity effect 
on the electronlike band. 

PACS numbers: 74.25.Gz, 74.20.Mn, 74.20.Rp 78.30.-j 



I. INTRODUCTION 

Pairing symmetry of electron-doped high-T c 
cuprate superconductors such as Nd 2 - 2: Ce 2 ;Cu04 
and Pr2- a; Ce a ;Cu04 is a long standing problem 

q Hi i ii ii Hi ii nam oi hi in Sim na. 

Although no consensus has been reached yet, more 
and more recent experimental results have sug- 
gested that the order parameter of electron-doped 
cuprates is li kely to have d x 2_ y 2-w&ve pairing symmetry 
[HM IS WiWi lid! IhH ITri Il7j . in close resemblance to that 
of hole-doped materials. Among various experiments, 
the Raman scattering seems to suggest a different story 
Q [H. For hole-doped superconductors, it was 
known that typical B\ g , Bi g , and A\ g pair-breaking 
peaks appear, respectively, at the frequencies of 2, 1.6, 
and 1.2 times of the gap amplitude 18]. However, in 
electron-doped materials the relative position of the B\ g 
and Bi g peaks changes with doping. The Big peak 
appears first at a higher frequency than the B\ g one in 
the underdoped regime. It then moves down and finally 
appears at a frequency lower than that of the B\ g peak 
in the heavily overdoped regime. 

The Raman scattering has the potential to probe dif- 
ferent regions of the Fermi surface (FS), thus a thorough 
understanding of the experimental data can provide a 
better understanding on the momentum dependence of 
superconducting (SC) pairing gap. The observation of 
Big Raman peak at higher frequency than that of B\ g 
would imply a non-monotonic ^2 ,,2 -wave order param- 
eter in a single-band system |l4| . This nonmonotonic 
order parameter seems to be also consistent with the ob- 
servation of angle resolved photoemission spectroscopy 
(ARPES) [3. However, this one-band picture may not 
be adequate to describe the nature of two kinds of charge 
carriers in electron-doped cuprate superconductors, as re- 
vealed by magneto-transport measurements [1!1 l2Cl l21| . 

A key clue towards the understanding of Raman data 
in Nd2- 3; Ce 3 ;Cu04 comes from the doping evolution of 



the FS revealed by ARPES |22j,|2J. At low doping, four 
small FS pockets first appear around (±7r, 0) and (0, ±n). 
By increasing the doping, four new pockets begin to form 
around (±7r/2, ±7t/2). These results can be explained in 
terms of the k-dependent band-folding effect due to anti- 
ferromangetic (AF) ordering. The original band is folded 
back into the magnetic Brillouin zone (MBZ) around the 
diagonal line (it, 0) — > (0, it). Near the intersecting points 
of the Fermi surf ace, an AF gap opens and splits the orig- 
inal FS into two |24l I25L |26| . This two-band picture was 
first used by Luo and Xiang to explain the unusual tem- 
perature dependence of the magnetic penetration depth 
in electron-doped copper oxides [23. It is supported 
by Hall coefficient and magneto-resistance measurements 
|19t l2Ci l21j . The generic feature of a weakly coupled two- 
band model was discussed in Ref. [28| in the context of 
hole-doped cuprate superconductors. 

In this paper, we shall use a two-band model to study 
the Raman response for the electron-doped cuprates. As 
will be shown, the two-band model gives a unified expla- 
nation to the unusual behaviors of Raman spectra. It 
explains in a natural way why the Bi g Raman peak ap- 
pears at a higher frequency than that of the B\ g peak at 
optimal doping and why the relative positions of these 
two peaks change in the heavily overdoped regime. 



II. MODEL AND FORMALISM 

We start by considering the two-dimensional t-t'-t"-J 
model 

H = -* £ 4^* - *' S 4,<*r - *" £ c *W 

<»i>2CT (»3>3CT 
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where (27)2, and (1,7)3 denote the nearest, second- 

nearest, and third-nearest neighbors between i and j. No 
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double occupied sites are allowed in Eq. (1). All nota- 
tions used in are standard. Applying the slave-boson 
and MF decoupling [2j|, Hamiltonian can be written 
in terms of two (diagonalized) bands in momentum space 

H = Y1 (^a^aica + £k/3^lLxAc<x), 

where the prime denotes that the momentum summation 
is over the MBZ only (— tt < k x ± k y < tt) and 

£ £k + gk+Q _^ / (£k+Q - £k) 2 , „ T2 2 

6c,i = ^ T y ^ + 4:J 2 m 2 

with Q = (tt, tt) the AF wave vector and m = (— 1) 1 {S' ) 
the AF order. Here 

gk = (2\t\5 — J x) (cos k x + cos k y ) 

— At' 5 cos k x cos k v — 2t"S(cos 2k x + cos 2k y ) (2) 

with x = (fiafjc) ^ ne uniform bond order and 5 the dop- 
ing concentration (fi a is the fermionic spinon operator). 

In the SC state, we add a BCS coupling term to each 
band and assume the system to be described by the fol- 
lowing Hamiltonian [27| 

H = ]T'^ k A- + H' Ak ^k T ^ki + i-k;Jkt), 

kal kZ 

where I = a, (3 and Ak,; = (A;/2) [cosfc^ — cosk y ] are the 
c? K 2_j / 2-wave gap functions. 

The Raman scattering intensity is proportional to the 
imaginary part of the effective density-density correlation 
function x(q, t) = (T T [p(q, t), p(— q, 0)]} in the limit q — ► 

0. Here p(q,r) = Ek,<T7kC k + q ( r ) c k( r ) is the effective 
density operator and 7k is the Raman vertex. When 
the energy of incident light is smaller than the optical 
band gap, the contribution from the resonance channel 
is negligible. The Raman vertex can then be obtained in 
terms of the curvature of the band dispersion under the 
inverse effective mass approximation. 

In the current two-band model, the effective density 
operator is decomposed as 

k,<r 

+ 7k+QC k +Q+q, CT ( T ) C k:+Q,<T(T)]. (3) 

Along with the unitary transformation |25| such that op- 
erators Ck and Ck+Q are transformed into ak and /3k+Q, 
the Raman response function for each symmetric channel 
(S) is then given by [2^ 

X S (q - 0,r) = -^'( 7 ^,) 2 [^(k,r)^(k,-r) 
k,W 

-c w ^,(k,r)^(k,-T)], (4) 



where Q and T are the normal and anomalous Green 
functions for a superconductor, ew = 1 if I = I' or — 1 if 
I I' . The intra- and interband vertex functions are 

7k,aa = cos 2 6^ +Sin 2 6>k7k+Q, 

7k,/3^ = sin 2 &k7k + c °s 2 ^kTk+Q) (5) 
7k,a/3 = 7k,/5a = sir^flk^ - 7k+q), 

where cos26> k = (gk+Q - gk)/V( £ k+Q - gk) 2 + 4J 2 m 2 . 
For the B\ g and B 2g channels, 7 k ?ls = — 7™ and 

7k 29 = ^J^- Here 7k J = d 2 e\ s _/dk i dkj (inverse effec- 
tive mass approximation). Eq. Q and JSJ reduce to the 
famous ones in a one-band system when m = Il8l . 
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FIG. 1: Schematic plot of the Fermi surfaces for the a and /3 
bands and the momentum dependence of the Bi g and Bi g Ra- 
man vertices. The B\ g and B29 modes couple more strongly 
with the a and /3 bands, respectively. 

Shown in Fig. ^ ar e the £?2g and Raman vertices 
and how they are coupled to the SC quasiparticle excita- 
tions in k space. Since the Bi g vertex has d xy symmetry, 
the B2 g channel is dominated by the excitations of the (3 
band. In contrast, the B lg vertex has d x 2_ y 2 symmetry, 
thus the Big channel is contributed mainly from the ex- 
citations of the a band. Since both B 2g and B\ g vertices 
are odd-parity, their Raman intensities are not affected 
by the Coulomb screening. For the fully symmetric A\ g 
channel, in contrast, all regions of momentum space con- 
tribute and the Raman intensity is partially screened. 
Since the A\ g channel is more sensitive to the actual ver- 
tex as well as the screening effect, we will leave out A\ g 
and focus on the B 2g and B\ g channels only. 

In electron doped cuprates, the SC state appears only 
when the [3 band emerges above the Fermi energy. This 
would suggest that we assume that it is the f3 band that 
drives the system to superconduct, while the a band be- 
comes superconducting mainly via the proximity effect. 
A simple picture for the understanding of the unusual 
Raman spectra in electron doped cuprates can then be 
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sketched as follows. In the under doped or optimally 
doped regime, the (5 band couples more strongly with 
the AF fluctuations than the a band. This results in a 
relatively larger SC gap in the (3 band (A/3) than in the 
a band (A a ). The B\ g channel probes mainly the quasi- 
particle (QP) excitations of the a band, thus the B\ g 
Raman peak is mainly determined by A Q . Similarly, the 
B 2g channel probes mainly the QP excitations of the j3 
band, its Raman peak is mainly determined by Ap. If A a 
is much smaller than Aa, one would then expect that the 
Big peak appears at a frequency lower than that of the 
Big peak, unlike in the hole doped case. In the heavily 
overdoped regime, the AF correlation becomes very weak 
and the band splitting vanishes. In this case the two- 
band model reduces essentially to a one-band model and 
the Big and Bi g Raman spectra would behave similarly 
as in the hole doped cuprate superconductors, consistent 
with the experiments. 



TABLE I: 


Summary of 


fitting 


parameters for 


Nd 2 -.C e;c Cu04 


(refer to Fig. 2) 








Column 1 


Column 2 Column 3 


T c (K) 


22 


22 


13 


m 


0.178 


0.178 


0.15 


T^cm" 1 ) 


6 


4 
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r^cm- 1 ) 


8 


6 


8 




21 


27 


21 


A^cm" 1 ) 


48 


57 


37 


Af,/A a 


2.29 


2.11 


1.76 


Ap/T c 


2.18 


2.59 


2.84 



high (low) energy region to which the Big (Bi g ) probes. 
It is thus expected that the Bi g {Big) peak will be red 
shifted (blue shifted) compared with the result of hole- 
doped cuprate superconductors. 



III. RESULTS AND DISCUSSIONS 

A. On Nda-^Ce^CuOi 

Pertaining to Nd2- K Ce a; Cu04, we have adopted |t| = 
0.326 eV, t' = 0M, t" = -0.2t, and J = 0.3t to sim- 
ulate the band structure j^. We take \ = —0.15 and 
m = 0.178 for the optimally doped (x — 0.15) sample, 
X = —0.15 and m = 0.15 for the overdoped (x = 0.16) 
sample. These parameters are close to those obtained 
in self-consistent calculations for the normal states [25| . 
The chemical potentials arc determined by the filling fac- 
tor for each band to give the true doping concentration 
through x = n e — rih- 

Theoretical fitting procedures are implemented as fol- 
lows. First, the vertex functions are evaluated using JSJ. 
Then the SC gaps A a and A 3 and the smearing Lorentz 
width r Q and Tp are adjusted to fit the peak positions 
and the overall spectral line shape (up to a constant mul- 
tiplying factor). 

Figure [21 compares the experimental data of Raman 
spectra for Nd2- a; Ce a: Cu04 with the theoretical calcula- 
tions. For all the cases considered in Fig. |3 our fitting 
curves are in good agreement with the experimental re- 
sults. For the optimally doped sample (x=0.15) reported 
in Ref. [3 at T = 8 K (first column), the B lg and B 2g 
Raman peaks appear at 50 cm -1 and 55 cm -1 , respec- 
tively. The corresponding gap and smearing parameters 
obtained by fitting are (A a ,Ap) = (21 cm -1 , 48 cm -1 ) 
and (T a ,Tp) — (6 cm" 1 , 8 cm" 1 ). The ratio between the 
Big Raman peak frequency and A a is about 2.4, while 
the ratio between the Bi g Raman peak frequency and Ap 
is about 1.2. The corresponding ratios in a hole doped 
d a; 2_j ( 2-wave superconductor are about 2 and 1.6, respec- 
tively. This difference between hole and electron doped 
cuprates is not difficult to be understood. In electron 
doped materials, the AF correlation splits the continu- 
ous FS into two separate sheets. This then suppresses the 
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FIG. 2: Comparison of the theoretical calculations 
(solid lines) with the measurement data (circles) for 
Nd2-a;Cea;Cu04. The experimental data in the first column 
are taken from Ref. |T3 |. The data in the second and third 
columns are taken from Ref.|l5jl. u)l is the incident photon 
energy. 

For another set of data of the optimally doped sample 
reported in Ref. [l5j at T = 4k (second column) , the fit- 
ting parameters are (A a , Ap) — (27 cm -1 , 57 cm -1 ) and 
(r Q , r^) = (4 cm , 6 cm" 1 ). The corresponding param- 
eters for the overdoped sample at T = 4 K (third column) 
are (A a ,Ap) = (21 cm" 1 , 37 cm" 1 ) and r Q = Yp = 8 
cm" 1 . The gap parameters obtained are consistent with 
the general expectation. Both A a and Ap decrease with 
increasing temperature at the same doping level and with 
increasing doping at the same temperature. The gap ra- 
tio, r = Ap/A a , is reduced from 2 at optimal doping to 
1.7 at slightly overdoping, consistent with the scenario 
of AF-like fluctuation induced superconductivity. All fit- 
ting parameters for Nd2_ a; Ce 2 ;Cu04 are summarized in 
Tabled 
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TABLE II: 


Summary 


of 


fitting 


parameters for 


Pr 2 - a! Ce a; Cu04 


(refer to Fi| 


3). 








Column 1 




Column 2 


Column 3 


X 


0.15 




0.165 


0.18 


T c (K) 


23.5 




15 


10 
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0.15 




0.12 
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8 


15 
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15 


AoK 1 ) 


31 




16 


15 


A^cm" 1 ) 


68 




30 


15 


A /A a 


2.19 




1.88 


1 




2.89 




2 


1.5 



B. On Prj-.CejCuOi 

The Raman scattering measurement has also been 
carried out in electron-doped P^-xCe^CuC^ at var- 
ious doping levels The Raman spectra of 
Pr2-a;Ce x Cu04, as shown in Fig. El behave similarly as 
for Nd 2 - x Ce x Cu04. At optimal doping, the B 2g peak 
appears at a frequency higher than that of the B\ g peak. 
With increasing doping, the frequency of the B\ g peak 
approaches to and finally surpasses the B 2g peak in the 
overdoped regime. 
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FIG. 3: Comparison between the theoretical results (solid 
curves) and the Raman measurement data from Ref. El (cir- 
cles) for Pr 2 _ 2 ,Ce :E Cu04 at T = 4K. lo l = 1.9eV. 

Figure. |21 compares the measurement data of Pr 2 _a; 
Ce^CuC^ with our theoretical calculations. For the op- 
timally doped sample (x = 0.15), the Raman peak ap- 
pears at 80 cm -1 for the B 2g mode and at 62 cm -1 for 
the B\g mode. The parameters obtained by fitting are 
(A Q , Ap) = (31 cm" 1 , 68 cm" 1 ) and T a — Tp — 6 cm" 1 , 
with x = —0.15 and m — 0.15. For the slightly over- 
doped sample [x = 0.165), the B 2g and B\ g peaks ap- 



pear at the same frequency at 37 cm 1 and the parame- 
ters we obtained are (A a , A^) = (16 cm" 1 , 30 cm -1 ) and 
r Q = = 8 cm" 1 , with \ = -0-15 and m = 0.12. For 
the heavily over-doped sample (x = 0.18), the Raman 
peaks appear at 25 cm -1 and 30 cm -1 for the B 2g and 
Bi g modes, respectively. The relative peak positions of 
these two modes are similar as in a one-band d 2 .2_ y 2-wave 
superconductor. This is not unexpected since at such a 
high doping level, the two-band model reduces essentially 
to a one-band model. In this case, the parameters we ob- 
tained are A Q = A^ = 15 cm -1 and T a = Tp = 15 
cm" 1 , with x = —0.15 and m = 0. For the above three 
samples, the ratio A^/A Q changes from 2.3, to 1.9, and 
finally to 1 with increasing doping. With increasing dop- 
ing, the AF order is depressed and the gap amplitudes 
is decreased. The results are consistent with neutron 
scattering measurement [3p| . All fitting parameters for 
P^-zCezCuC^ are summarized in TABLE ITU 

In the top-left panel of Fig. [3] for the B 2g mode, the 
experimental data are taken under the strong resonant 
regime as emphasized in Refs. flil Il5j . In this case, 
the contribution from the resonance channel becomes im- 
portant and the inverse effective mass approximation is 
not valid. Our theoretical result including only the non- 
resonant contribution can give a good explanation to the 
low-frequency part of the spectrum, but the height of the 
peak is much lower than the experimental one. The res- 
onance channel may also have some contribution to the 
B 2g spectrum in the second column of Fig. [2] for Nd2-a; 
Cea;Cu04. 



IV. SUMMARY 

In summary, we have analyzed the Raman spectra 
of electron-doped cuprate superconductors based on a 
weakly coupled two-band model. Our result gives a uni- 
fied explanation to the experimental data in the whole 
doping range. It suggests strongly that the SC pairing in 
electron-doped cuprate superconductors results from the 
same pairing mechanism as in hole doped ones and the 
gap parameter has d a .2_ y 2-wave symmetry. To under- 
stand the Raman data in the strong resonance regime, 
a more comprehensive theory including the contribution 
from the resonance channel is desired. 
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